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Abstract 

In this paper, we give a muhistep extension of the e-algorithm of Wynn, and we show 
that it implements a multistep extension of the Shanks' sequence transformation which is 
defined by ratios of determinants. Reciprocally, the quantities defined in this transformation 
can be recursively computed by the multistep e-algorithm. The multistep e-algorithm and 
the multistep Shanks' transformation are related to an extended discrete Lotka-Volterra 
system. These results are obtained by using the Hirota's bilinear method, a procedure quite 
useful in the solution of nonlinear partial differential and difference equations. 



1 The scenery 



Let {Sn) be a sequence of numbers converging to S. If its convergence is slow, it can be trans- 
formed, by a sequence transformation, into a set of new sequences {(T"^"^)}, depending on two 
indexes k and n, and converging, under certain assumptions, faster to the same limit, that is such 
that 

rjnin) _ g rp{n) _ g 

lim — = 0, or lim — = 0, or both. 

A well-known example of such a transformation is the Richardson extrapolation process, 
which gives rise to the Romberg's method for accelerating the convergence of the trapezoidal rule 
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for approximating a definite integral. Let us mention that sequence transformations can also be 
applied to diverging power sequences, thus leading, in some situations, to interesting results such 
as analytic continuation (this is the case of the e-algorithm which, applied to the partial sum of 
a divergent power series, computes its Pade approximants). 

In many sequence transformations, the terms of the new sequences can be expressed as ratios 
of determinants, and there exists, in each particular case, a (usually nonlinear) recursive algorithm 
for avoiding the computation of these determinants and implementing the transformation under 
consideration [IIl[38lll2Hll] • 

The most well-known transformation of this type is due to Shanks [361137] . It can be imple- 
mented via the e-algorithm of Wynn [l5]. Recently, a new recursive algorithm for accelerating 
the convergence of sequences was derived by He, Hu, Sun and Weniger [T3] from the lattice 
Boussinesq equation. This algorithm resembles to the e-algorithm, and it was proved that the 
quantities it computes can be expressed as ratios of determinants, thus extending the Shanks' 
sequence transformation. In this paper, inspired by this approach, we will extend further the e- 
algorithm, and we will show that it implements an extension of the Shanks' transformation, thus 
leading to a multistep e-algorithm and a multistep Shanks' transformation. The proof makes use 
of the Hirota's bilinear method [T7] which was invented for resolving integrable nonlinear partial 
differential or difference evolution equations having soliton solutions. 

For some years now, there has been a great concern for convergence acceleration algorithms 
among the community of mathematical physicists working on integrable systems, KdV and other 
equations, soliton theory, Toda lattices, etc. [9l|2ll[25l[3T],|32]. These researchers are interested 
by the fact that convergence acceleration algorithms are nonlinear difference equations in two 
variables whose solutions are explicitly known. Determinants often play a central role in this 
type of problems as exemplified, for example, in [H]. An important procedure for obtaining a 
closed-form solution of soliton equations is the Hirota's bilinear method [T7] which consists in 
writing the solution as a ratio, and then working with its numerator and its denominator. 

In Section [21 we discuss the Shanks' sequence transformation and its implementation by the 
e-algorithm of Wynn. The quantities involved in this transformation and in this algorithm are 
expressed by ratios of Hankel determinants. In Section [3l we present our multistep extension of the 
^-algorithm, and the corresponding multistep extension of the Shanks' transformation. Section H] 
is devoted to some relations between determinants that will be useful for our purpose. The Hirota's 
bilinear method is presented in Section O In Section El we show that the quantities recursively 
computed by the multistep e-algorithm correspond to the ratios of determinants defining the 
multistep Shanks' transformation, and, reciprocally, in Section [71 we show that the multistep 
Shanks' transformation can be implemented by the multistep e-algorithm. Finally, in Section [HI 
the connection between an extended discrete hungry Lotka-Volterra system and the multistep 
e-algorithm is discussed. Hirota's method is essential for obtaining these results. The paper ends 
by some considerations on further researches. 
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2 The Shanks' transformation and the e— algorithm 



The Shanks' sequence transformation [36l[37] : (Sn) i — > {{^kiSn))} consists in transforming a 
given sequence (Sn) into the set of sequences {(efc(S'„))} whose terms are defined by 



ek{Sn) 



k,n = 0, 1, 



where A is the usual forward difference operator whose powers are defined by 
with A^Sn = Sn, and where Hkiun) denotes the Hankel determinant 



(1) 



Ur. 



Un+l 



Un+k~l 



with 7io{un) = 1- 

Obviously, replacing each row, in this determinant, by its difference with the previous one, 
repeating this operation several times, and performing it also on the columns, we have 



T-LkiUn) 



Un 
AUn 
A^Ur, 



A'^^^U 



Un+k-1 
AUn+k-1 

■n+k-1 



A^u 



■ A^ ^M„+fc-l 



Un AUn 
AUn A^Un 



A^Un 



2k-2„ 



The e-algorithm is a recursive algorithm due to Wynn [U] for implementing the Shanks' 
transformation without computing the Hankel determinants appearing in ([1]). Its rule is 



(n+1) 



^k+1 — ^k-l 



(n+1) (n) 



k,n = 0, 1, 



e 



(2) 



el 



with e^^l = and e^^^^ = Sn, for n = 0, 1, . . . 

The connection between the e-algorithm and the Shanks' transformation is given by 



^2k — ^k{Sn) and £:2fc+i 



(n) 



1 



ek{ASn 



k,n = 0, 1, 



(3) 



Thus, the £^2^+1'^ intermediate results, and we have 

(n) _ 'Hk+ljSn) , (n) _ TikiA^Sn) 

- HkiA^Sn) ~ m+liASn) 



The quantities e]^ are usually displayed in a two-dimensional array (the e -array) where the 
lower index k remains the same in a column of the table, and the upper index n is the same in 
a descending diagonal. Thus, the rule (|2]) relates four quantities located at the four vertices of a 
lozenge in three different columns and two descending diagonals as showed below 

k 

(n+1) (n) 
(n+1) 

For implementing the e-algorithm efficiently, the best technique, due to Wynn [ITjHS] . consists 
in storing the last ascending diagonal of the e-array (in this diagonal the sum of the lower and the 
upper indexes is constant), and to add, one by one, the terms of the sequence to be transformed. 
Then, a new ascending diagonal is built step-by-step, by moving up the lozenge, and the new 
diagonal gradually replaces the old one. The corresponding FORTRAN subroutine can be found 
in [n]. 

Since the quantities with an odd lower index are intermediate computations, they can be 
eliminated, thus leading to the cross rule also due to Wynn [12] 

1 1 _ 1 1 

~{n) _ ^(n+1) + (n+2) _ (n+1) " (n+2) _ (n+1) + ~H _ Jn+l) ' 
^2fc+4 ^2fc+2 ^2k ^2fc+2 ^2fc+2 ^2/fc+2 ^2fc+2 ^2fc+2 

with the initial conditions = ^ ^^'^ ^o""* = 5'„ for n = 0, 1, . . . Obviously, it is also possible to 
eliminate the ^["^'s with an even lower index for obtaining a rule only involving quantities with a 
lower odd index, although this is less useful from the numerical point of view. 

The proof given by Wynn for his e-algorithm was mostly a verification of the link between the 
Shanks' transformation and the algorithm, since he introduced the ratios of Hankel determinants 
for ek{Sn) and efc(AS'„) into the rule of the e-algorithm, and he showed that the equahty held by 
making use of the Sylvester's determinantal identity and the Schweins' one which can be found, 
for example, in (see [SI pp. 142-143] for their proofs). The difficulty of the proof resided in 
the nonlinearity of the algorithm. Of course, Wynn's great merit was the idea of the e-algorithm 
itself, followed by this verification. 

There are three approaches for linking a sequence transformation and a (usually nonlinear) 
recursive algorithm for its implementation. By increasing order of complexity, they are 

1. Verification: the transformation and the algorithm are both known, and one has to verify 
that they lead to identical sequences. This is the way followed by Wynn in [15] when he 
gave his e-algorithm. 

2. Derivation: only the transformation is known, and one has to derive an algorithm for 
its implementation. This is the case, for example, of the ^'-transformation which is the 
most general sequence transformation known so far, and which can be implemented by 
the ii^-algorithm, an algorithm which appeared almost simultaneously in various contexts 
[HI [131 [221 135] . This was also certainly the way Wynn followed when he derived his e- 
algorithm, although it was not presented like that in his paper |45] . 
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3. Proof: only the algorithm is known, and one has to guess a formula (that is a ratio of 
determinants) for the transformation it is implementing, and to prove it. This was the 
situation for the second generalization of the e-algorithm proposed in [5], whose form was 
obtained by Salam [33|IM]. Let us mention that the ^-algorithm [Ij is an extrapolation 
algorithm for which no determinantal formula is known yet, if it exists. 

Now, after presenting the multistep e-algorithm and the multistep Shanks' transformation 
(Section [3]), we will show, with the help of determinantal identities (Section Hj) and the Hirota's 
bilinear method (Section |5]), how to go from the multistep e-algorithm to the multistep Shanks' 
transformation (Section [6]), and back (Section [7]). 

3 The multistep £— algorithm and the multistep Shanks' 
transformation 

Let m be a fixed strictly positive integer. We define the multistep e -algorithm by the recursive 
rule 

(n) _ (n+1) 1 L Q 1 (A\ 

with the initial values 

_ n - - ... - - „ r^") - S* « - 1 f^') 

Displaying these quantities in a double array similar to the e-array, we see that this rule relates 
2m + 2 quantities located in an extended lozenge covering m + 2 columns and two descending 
diagonals as showed below 

An) 
k—m+l,m 
(n+1) (n) 
k—m,m k—m+2,m 

fc— m+l,m 

M 

^k,m 

Jn+l) in) 

'^k—l,m ^k+l,m 
(n+1) 



^k,m 



The implementation of this algorithm using the technique of ascending diagonals, as described 
above for the e-algorithm of Wynn, is more difficult, and it requires the storage of m ascending 
diagonals for computing the (m + l)th one. 
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In Section [6], we will prove that, for all k and n, it holds 

in) 



'(m+l)(fc-l)+l,m 



, 2 = 2, 3, ...,m, 



-(m+l)(fe-l)+i,m HkiA'Snj 

where the determinants Hk and which depend on m, are defined by 



(6) 
(7) 
(8) 



HkiUn) 



Un 
^2m 



n+l 



Ur. 



A 



n+l 



2m 



U 



n+l 



A'^Un+k-l 
A'^"^Un+k-l 



with H-iiun) = and Ho{un) = 1, and where 



1,2, 



n = 0, 1, . . . , 



A-u„ 



^2m 



n + l 
A'^u^+i 

Un+1 



^2m 



n + k — 1 

A^Un+k~l 
Un+k-1 



^2m 



^{k-2)r. 



^{k-2)m^ 



n+l 



^(k-2)r. 



Un+k-l 



1,2, 



n 



0,1,... 



with <l>_i(n„) = and ^o{un) = 1- 

Let us notice that, when m = 1, Hk{un) is identical to the usual Hankel determinant Hkiun)- 

For proving these determinantal identities, we will follow a procedure similar, although more 
difficult, to the procedure used in [T?j (which is based on the Hirota's bilinear method) for deriving 
a determinantal expression for a new acceleration algorithm obtained from the lattice Boussinesq 
equation. However, instead of the Jacobi's determinantal identity, we will only use the Sylvester's 
one (which is, in fact, the same after a permutation of rows and columns), and we will not use 
the Schwein's identity. 

Let us now define the multistep Shanks' transformation ek,m '■ (Sn) ' — > {{Gk,m{Sn))} by 



,m 



-(«) 



{m+l)k,m HkiA'^+^Sn 

Obviously, after proving ([7]), we will also have 



/c, n = 0, 1, 



(9) 



'(m+l)fc+l,m 



1 



,m 



fc, n = 0, 1, 
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a result similar to the second relation ([3]) for the e-algorithm of Wynn. Thus, only the quantities 
^fcm'^ whose first lower index is a multiple of m + 1 are interesting for the purpose of convergence 
acceleration. All the other ones are intermediate computations. The computation of Ck^mi'Sn) = 

^(m+i)fc,m ^^eds the knowledge of S'n, • • • , Sn+mk- 

For simplicity, we will omit to indicate that all the symbols used in this paper depend on the 
fixed integer m. 

We see that, when m = 1, the algorithm (j4]) reduces to the ^-algorithm (|2]), and the trans- 
formation (Q reduces to the Shanks' transformation ([1]). When m = 2, the recursive rule (j4]) 
reduces to the algorithm obtained in [H] from the lattice Boussinesq equation; see also [2^130] . 

Let us mention that, due to ([6]), the multistep Shanks' transformation can likewise be imple- 
mented by the E'-algorithm [B] with gi{n) = A*™S'„ for i = 1,2,..., and for all n, and that we 
get, for all k and n, E^^ = ek,m{Sn)- Thus, by the fundamental property of the i?-algorithm, 
the kernel of the transformation (jH]) (that is the set of sequences which are transformed into a 
constant sequence) is given by the 

Theorem 1 

A necessary and sufficient condition that, for all n, Ck^miSn) = S is that there exist constants 
Oi, . . . , Ofc, afc 7^ 0, such that, for all n, 

Sn = S + aiA'^Sn + asA^-^, + ■ ■ ■ + afcA'"™^„. 

Let us remind that the kernel of the Shanks' transformation Ckm '■ (Sn) ' — > {ekm{Sn) = ^2^™) 
is the set of sequences such that, for all n, Sn = S + biASn + ■ ■ ■ + hkmA^'^Sn, where hi, ... , bkm, 
bkm 7^ 0, are constants. Thus, we have the 

Corollary 1 

The kernel to the multistep Shanks' transformation e^.m is contained into the kernel of the Shanks' 
transformation Ckm- 

Moreover, due to the connection with the E'-algorithm, all the convergence and acceleration 
results proved for it OET] also hold for the multistep Shanks' transformation. 

In the next Sections, we will link the multistep Shanks' transformation ([9]) and the multistep 
^-algorithm (j3j) by means of the Hirota's bilinear method. First, in Section |U some relations 
between the determinants Hk{A^Sn) and <I>fc(A*S'„) will be established. We will only employ 
the Sylsvester's determinantal identity, contrarily to the proofs given in [IS] and [Hj where the 
Schweins' determinantal identity is also used. Then, Hirota's bilinear method will be presented in 
Section [5l In Section [6l we will show that the quantities computed by the multistep e-algorithm 
(jlj) are those defined in the multistep Shanks' transformation ([6])-([8j). Conversely, in Section 
[TJ we will prove that the multistep Shanks' transformation ©-([H]) can be implemented by the 
recursive rule (jl]) of the multistep e-algorithm. 
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4 Relations between determinants 



Let A be a square matrix, a,/?, 7 and 6 numbers, a,b,c and d vectors of the same dimension as 
A. Let M be the matrix 

(a a" 
b A c 
7 c/^ S 

The Sylvester's determinantal identity is 



T 

a a 




A 


c 




(3 




b 


A 


b A 






5 




A c 




7 


d^ 



\M\ ■ \A\ 

Let us now prove some determinantal identities that will be useful in the sequel. 
Lemma 1 

Proof: we consider the determinant 



(10) 



n+1 



A'S, 



n 

S„. A'+'^^Sn+l 



A^Sn+1 



n+k+l 



Ai+km Q 
Jn+k+1 

A^Sn+k+l 



{-irHk+^iA^^'S, 



The second expression for Di is obtained by replacing each column, from the last one, by its 
difference with the previous one. Thus, we get a determinant whose first row only contains except 
in the first column where the element is equal to 1. Expanding this determinant with respect 
to its first row, and putting its last row as the first one, we see that Di = {—l)''Hk+i{A^'^^Sn)- 
Let us now apply the Sylvester's identity to the first expression of Di, and perform a similar 
manipulation on the rows and the columns of the other determinants, we obtain 



n+l) 



iJ,(A^+-+i^„+i)i7,+i(A*^0. 



Setting z = in this relation, we get ffTUl) . ■ 

A similar identity, which will be used in the sequel, also holds if Sn is replaced by AS„ 

Lemma 2 



Hk^liA'^'Sn)HkiA'Sn+l) - Hk-liA'^'Sn+l)HkiA'Sn). (11) 

Proof: let D2 be the determinant obtained from Di by replacing k hj k — 1, and moving the last 
row to the second position. Replacing each column, from the last one, by its difference with the 
previous one, we see that D2 = iJfc(A*+^S'„), and, applying the Sylvester's identity to it, we get 
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(12) 



Lemma 3 

Hk{/S.Sn)Hk{/y" Sn+l) = Hk{/V^~^^ Sn)Hk{Sn+l) — -f^fc+l (5'„)iJfc_i ( A"""*"^ ^n+i ) . 

Proof: Setting z = m in flTT]) . we have 
Hk{A"'+'Sn)Hk^iiA^Sn+i) = Hk-i{A^^^Sn)Hk{A^Sn+i) - i/fc_i(A'"+^S„+i)i/fc(A'"S„). (13) 
Applying now the Sylvester's identity to the determinant i/fc+i(A*S'„,), we get 

-f^A:+l(A*S'„,)ifjt_i(A*^'"S'„+i) = Hk{A^ Sn)Hk{A^~^'^ Sn+i) — Hk{A^ Sn+l)Hk{A^^"^ Sn) ■ 

Setting i = in this relation, we obtain 

Hk+\{Sn)Hk-i{A'^ Sn+i) = Hk{Sn)Hk{A"^Sn+i) — H k{S n+l) H k{A"^ S n) ■ (14) 

Then, we multiply (IT^ by Hk{Sn+i), we multiply (fn|) by Hk^i{A"^^^ Sn+i) , and we subtract. 
It gives 

Hk^liA^Sn+l)[Hk{A^~^^Sn)Hk{Sn+l) — Hk+l{Sn)Hk^l{A'^~^'^ Sn+l)] 
= Hk{A^ Sn+l)[Hk-l{A'^~^'^ Sn)Hk{Sn+l) — Hk{Sn)Hk-l{A^~^^Sn+l)]■ 
\Jsmg (ITOl) . we see that the bracket in the right hand side is equal to Hk{ASn)Hk-i{A"^ Sn+i) ■ 
After simplifying both sides by Hk^i{A"^Sn+i), we obtain (IT^ . I 

A similar identity, which will be used in the sequel, also holds if 5"^ is replaced by ASn- 
Lemma 4 



HkiA'Sn+l)Hk^liA'^'Sn) = HkiA''-'Sn)<^kiA'Sn+l) - i/fc-l(A^+^^„+l)$fc+l(A*^0. (15) 



Proof: we consider the determinant 

1 



A^Sn 
A'+'^Sr, 



1 

A'S, 
A^+™5, 



n+l 

n+1 



n n + 1 



Obviously, we also have 



1 

n 
A'S„ 



1 

n + 1 
A^^„+i 



l)m^ /\i+{k~l)'m q 



n+1 



A'+^S„ 



1 

A^+15, 



1 

A^Sn+k+l 
A^^"^ Sn+k+1 

A i+(fc-l)m c 

n + k + 1 



1 

n + k + 1 

A^Sn+k+l 
1 



n+1 



^i+l+(k-l)mg^ ^i+l+(k-l)mg^^^^ 
-l)'Hk{A^^'S^). 



^ '~>n+k 
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We apply now the Sylvester's identity to the first expression of given above, and replace D3 by 
(— l)'^i/fc(A*"'"^5'„). We get, after similar manipulations on the columns of the other determinants. 

Then, we apply the Sylvester's identity to the determinant ^k+iiA^ Sn). We get 

<l>fc+i(A^S'„)iffc„i(A*5„+i) = MASn)HkiA'Sn+i) - MASn+i)HkiAS^). 

We multiply this identity by iJyfc_i(A*+^5'„+i), we multiply ffTTl) by $jt(A*S'„,+i), and we sub- 
tract. It gives 

Hk-l{ASn+l)[<^k+l{ASn)Hk^^{A+'Sn+l) " //fc ( A^+^ 5„) <l>fc ( A*5„+i)] 

= i/fc(A^S'„+i)[$fc(A^5„,)iffe_i(A^+i^„,+i) - i/fc_i(A^+i5„)$fc(A^5„+i)]. 

Using f|T6|) . we see that the bracket in the right hand side is equal to — iffe„i(A*"''^S'„)i7fc-i(A*5'„+i). 
After simplifying both sides by Hk^i{A^ Sn+i) , we obtain f|T5l) . I 

Lemma 5 

Hk{ASn)Hk^2{A''+'Sn+l) = Hk.^{A''+'Sn+l)Hk-l{ASn) - Hu^^{A^^'Sn)Hk-l{ASn+l). (17) 

Proof: we consider the determinant 



Sn 
1 



Sn+1 
1 

A^Sn+1 



Sn+k 
1 

A m c 
^ ^n+k 



After exchanging the first row and the second row, we see that D4 = —Hk{ASn)- Let us now 
apply the Sylvester's identity to D4, and perform a similar manipulation on the first and second 
row of the various determinants. We obtain ffTTl). I 



5 The Hirota's bilinear method 

The Hirota's bilinear method [T7] is a technique which could be much useful for solving certain 
nonlinear differential and difference equations. It consists in expressing the unknown as a ratio 
and, then, in treating separately the numerator and the denominator. 
We will now apply this method to the multistep e-algorithm, and set 



'k,m 



k 



We first have the 
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Lemma 6 

m m 

l-^fc+m+l'-^fc '^k+m+l) llK-'^k+i^k+i ^k+i^k+i)— fe+l-^fc+m 1 1 ' \^^) 

i=l 1=1 

Proof: 

Plugging f lTSjl into the recursive rule (jl]) of the e-algorithm, we get 

/^n f^n+l -■ 



k+1 k—m TT™ / k—m+i ^ k—m+i 

rpn+1 pn 
k—m+i k-m+i 



UZi ( 



-"^k-m^k+l -'^k+l^k-m _ I U=l k-m+i-'^ k-m+i 



pn pn+1 T-rm / p„, _ pn+1 r^n V 

^ k+l^ k-m lli=l\^ k-m+i^k-m+i ^ k-m+i^ k-m+il 



(20) 



Now, we cross-multiply the numerator of one side by the denominator of the other side, and 
we equate both sides. Replacing khy k + m and changing the sign, the equation fl2Up becomes 
(fig]) since 

m m+1 m 

Tpn I I TP'"- p^+l I I pn I I p»i+l ■ 

^k+m+l^k \_\_^k+i^k+i ~ 11 -^fc+j XX -^fc+i • ■ 

i=l i=l i=0 

The second preliminary result is contained in the 
Lemma 7 

// the following relations hold 

Tpn /^n+1 P"+l P" ('91 



(m+l)fc+l'-^(m+l)fc+l (m+l)fc+l (m+l)fc+2-' (m,+l)fc' 

pn /^n+1 pn+1 /^n rpn rpn+1 

{m+l)k+l^ (m+l){k-l)+l -^(m+l)(fc-l)+l'-^Cm+l)fe+l ~ (m+l)(fe-l)+2-^ {.m+l)fe ' 

and, for i = 2, . . . , m + 1, 



pn y^n+1 pn+1 /^n pn pn+1 ('9'?'l 

-^(m+l)fc+i'-J'(m+l')fc+j -^fm+l)fc+j'-^(m.+l)fc+i ~ (m+l)fc+j+l-^ (m+Dfe+i-l ' l'^'-'^' 



(m+l)fc+i'-^(m+l)fc+j {m+l)k+i^ {m+l)k+i ^ (m+l)fc+i+l-' (m+l)fc+i-l' 

pn /^n+1 pn+1 rin pn pn+1 

(m+l)k+i^ (m+l){k-l)+i ^ {m+l)(k-l)+i^ (m+l)k+i ~ {m+l){fc-l)+i+l (m+l)fc+i-l ' 

then f lT9|l follows. 



Proof: 

Let us first notice that, taking i = 1 in (12^ and gives (1^ and fl22]) . respectively, after a 
change in the signs of their right hand sides. 

Let us separate f lT9|) into the product of two relations, and prove that each of the following 
formulae holds separately 

m m 

U(^k+^Gltl - F-^lGlJ = ± n Fk+^+lFktU (25) 

i=l 1=1 
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which are the products appearing in both sides of (fT9l) . and 

which are its remaining parts. Then, multiplying together ( 125|) and (l26l) . we will obtain (fT9|) . but 
we must notice that the signs used in (125]) and have to be opposite. 

Let us assume that fl21l) - fl24p hold true. The proofs of the relations (1251) and (12^ have to be 
separated into three cases according to the value of k in (IT^ . 

• k replaced by (m + l)k in (125 p and (126 p . 

o Multiplying together the relations ( 123|) for z = 2, . . . , m, and then multiplying each of its 
sides by the corresponding side of (|2T]) (which brings a change in the sign) proves (|25|) . with the 
sign — , when /c is replaced by (m + l)/c in ( 125|) . 

o Replacing k hj (m + 1)A; in ( 12B]) . with the sign +, gives 

pn rpn+l /-in Tpn Tpn+1 

-^(m+l)(fc+l)'-^{m+l)fc ^ {m+l)k^ (m+l){k+l) " (m+l)fc+l-^ (m+l)fc+m' 

which is (IM|) for z = m + 1 . 

o We get (|T9|) by multiplying together the two relations. 

• k replaced by (m + l)k + 1 in ( 125|) and ( 126|) . 

o Multiplying together the relations ( 123|) for z = 2, . . . , m + 1. The result is the same as adding 
1 to all the lower indexes, and making the product for z = 1, . . . , m, which is (125]) with the sign 
+. 

o Replace by (m + + 1 in ( 126]) . with the sign — . It is exactly ( 122]) with k + 1 instead of 

k. 

o Multiplying together the two relations, we obtain (|T9]) . 

• k replaced by (m + l)k + j, for j = 2, . . . , m, in (125]) and (126]) . 

o Let 2 < j < m be fixed. Let us write that ( 123]) holds with 2<i+j<m + l instead of i, 
that is for z = 1, . . . , m — j + 1, 

Tpn /~in+l jpn+l /~in rpn rpn+l 

^ {ni+l)k+i+j^ {m.+l)k+i+j ^ {m+l)k+i+j^ {m+l)k+i+j ~ ^ (m+l)k+i+j+l^ {m+l)k+i+j-r 

Multiply together these relations for i = l,...,m— j + and, then, make their product 
for i = m — j + 2, . . . ,m. When i = m — j + 2, the left hand side of this expression be- 
comes Fl'^+i)k+m+2G'l;^ii)k+m+2 ' ^(^+i)fc+m+2^fm+i)fc+m+2 ' and its right hand side is equal to 
-^(m+i)fc+m+3-^{m+i)fc+m+i' "^^^t is, respectively, Fll^^-^^^(^f,^-^^^^G^^^^^^f^f^^^^^^—Ff^^^_^_^^f^f^^^^^ 
and -f'(m+i)(fc+i)+2-^(^+i)(fe+i)- Thus, by (|2T]) with k replaced by /c + 1, these two expressions are 
equal after changing the sign in one side. For z = m — j + 3, we have 



rpn /~in+l fpn+l /~in rpn jpn+l 

-^(m+l)fc+m+3'-^(m+l)fc+m+3 ^ {m+l)k+m+Z^ {m+l)k+m+'3, ~ (m+l)fc+m+4-^ (m+l)fc+m+2 ' 



that is 



-pn /^n+1 Tpn+1 /~tn pn Tpn+\ 

-^(m+l)(fc+l)+2'-^(m+l)(fc+l)+2 (m+l)(fe+l)+2'-^(m+l)(fc+l)+2 " -^^ (m,+l)(fc+l)+3-^ (m+l)(fc+l)+l ' 
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which is ( 123|) with k + 1 instead of k. And so on until i = m. Thus in the products from 
i = m— j + 2 to m, the sign is changed in one, and only one, of the expressions due to ([21]), and 
we finally obtain f l25p with the sign — . 

o Let us replace k hj {m + l)k + j in ( 126|) . with the sign +. We get 

^ {m+l)(k+l)+j^ {m+l)k+j ^ {m.+l)k+j^ {m+l)(k+l)+j — ^ {m+l)k+j-^ (m+l){k+l)+j~l^ 

which is when k is replaced by k + 1. 

o The product of the two relations gives fll9p . 

Thus, f ll9p have now been proved for all values of k. ■ 
Finally, we are able to prove the 

Theorem 2 

The relation f[T^ holds with the 's and the 's given by the following relations, for k = 
0,1,..., 

= H,{A'Sn), z = l,2,...,m + l (27) 

^Cm+l){A:-l)+l = Hk^i{A"^^'^Sn), C^^^.^^^ = iJfc+i(5'„), (28) 

= $fc+i(A-i5„), z = 2,3,...,m. (29) 

Proof: 

We are now able to prove ([2l])-(l24]), with the F^'s and G^'s given by ([27])-([29]). Replacing the 
determinants in ( ITOj) by their expressions, we obtain 

^ {m+l){k+l)^ {m+l){k+l) ^ {m+l){k+l)^ {m+l){k+l) " -^(m.+l)(fc+l)+l-^(m+l)fe+m 

which corresponds to fl25]) for the case z = m + 1. Replacing the determinants in flTB]) by their 
expressions, we obtain the bilinear equation ( 123|) for the cases i = 2,3, ... ,m, which completes 
the proof of the equation ( 123|) . 

Replacing the determinants in ( ITOj) and ( IT2|) . both with AS'„ instead of S'„, by their expressions, 
we see that the equations (ET]) and ([22]) are satisfied. 

Then, replacing the determinants in f ll2p by their expressions, we obtain 

^ {m+l)k^ {m+l){k~l) ^ {m+l){k-l)^ {m+l)k " -^^ (m+l)(fc-l)+m 

which corresponds to ( 124)) for the case z = m + 1, while replacing the determinants in (fT5]) by 
their expressions, we get the bilinear equation (121)) for i = 2, 3, . . . , m, which completes the proof 
for the equation (121)) . 

Since the identities ([2l])-([2l]) hold, then ([I9]) follows with the F^'s and the G^'s given by 
(ED)- ([29]). I 

We also have the 
Corollary 2 

-'^fc+m+l-'^fc-l ~^k^k+m ^k+m^k 
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Proof: 

Replacing the determinants in (ITT]) by their expressions given by ( HT^ . and k hj k + 1, we obtain, 
for i = 1, . . . , m, the following relation without any 

-^{m+l)fe+i+l-^(m+l)(fc-l)+i ~ -^(m+l)(fe-l)+j+l-^(m+l)A:+j ^ {m+l)k+i^ (m+l){k-l)+i+V 

Similarly, the determinantal identity f[T7|l leads, after replacing /c by A; + 2, to 

-^(m+l){fc+l)+l-^(m+l)fc ~ -^(m+l)fc+l-^(m+l)(fc+l) (m+l)(fc+l) (m+l)fc+l ' 

which is the preceding relation for i = m + 1. Thus, changing mk into m, these two identities 
can be gathered into the single formula of the Corollary. I 



6 From the multistep algorithm to the multistep Shanks' 
transformation 

By comparing (l6])-([8]) with the determinantal formulae (I27|) - (l29|) of the Theorem [2] for the -F^'s 
and the G^'s issued from the Hirota's method, we are now able to give the determinantal formulae 
for the multistep e-algorithm. Consequently, from the Lemmas El [71 and the Theorem [21 we have 
the 

Theorem 3 

The quantities e^*^^ computed by the multistep e-algorithm with the initializations (jSj), are 
expressed by the ratios of determinants (jSj); dZD, cin'd 

Thus, starting from the determinantal identities between Hk{A^Sn) and $fc(A*S'„), we proved 
that ([2ID-(I21D are satisfied with the determinantal formulae ([27D-([29D for the FJ^'s and the G^'s. 
Then, (|T9l) followed, and we concluded that the determinantal expressions ([6|)-([8|) for the e[."m's 
hold true. Notice that all these results were obtained without using the rule (|4]) of the multistep 
£-algorithm. 

Let us remind that, as noticed in [2] and fully explained in [7], we have 

^(n) _ fk,m{,Sn,---,Sn+{m+l)k) ^(n) _ D fk^m{ASn, ■ ■ ■ , ASn+{m+l)k) 

(m+l)fc,m Dfk^rn{Sn, ■ ■ ■ , Sn+{m+l)k) (m+l)k+l,m fi^ ,^(^ASn, • • • , ASn+(m+l)k) 

where fk,m is a function depending on (m + l)k + 1 variables and such that D'^fk^m = 0, where 
Dfk,m denotes the sum of the partial derivatives of fk,m- Thus, we obtain the following connection 
with Hirota's bilinear method 

G'(m+l)k ~ fk,m{Sn-, ■ ■ ■ i Sn+{m+l)k)i -^{m+l)fc ~ Dfk^m{Srn ■ ■ ■ i Sn+{m+l)k)-! 

^(*m+l)fc+l ^ Dfk,m{ASn, • • • , A5'„+(m+l)fc) ' -^(m+l)fc+l ^ fk,Tn{ASn, • • • , ASn+{m+l)k), 

and, according to this theory, the multistep Shanks' transformation is quasilinear that is ek,m{0'Sn+ 
h) = aCk^miSn) + b, a result which can be seen directly from ([9]). 
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7 From the multistep Shanks' transformation to the mul- 
tistep £— algorithm 

We will show now how to derive the recursive rule (jl]) of the multistep e-algorithm from the 
definition ©-([H]) of the multistep Shanks' transformation. 
From the determinantal identity f lT^ . we get 



An) _ (n+l) 

*^(m+l)(fc+l),m ^(m+l)fc,m 



Hk+2{Sn) Hk+l{Sn+l] 



Hk^iiA"-+^Sn)Hk{A--+'Sn+i) 

Hk+l{ASn)Hk+l{A'^Sn+l) 



i/fe+i(A-+iS„)iJ,.(A-+iS„+i)" 
Similarly, by the identity (IT^ with Sn replaced by ASn, we get 

in) ('^+1) 

fffc+i(A™+25„) Hk{A"^+'Sn+i: 



(30) 



'(„+l)(fc+l)+l,„, -^m+l)k+l,m Hk+2{ASn) Hk+l{AS, 



n+l) 

Hk+ljA'^'^'^ Sn)Hk+l{ASn+l) — HkjA^^'^'^ Sn+l)Hk+2{ASn 

Hk+2{ASn)Hk+l{ASn+l) 
Hk+l{A'^ Sn)Hk+l{A'^~^^ Sn+l) 



Hk+2{ASn)Hk+l{ASn+l] 

We also get the following relation from the identity ( IT5|) 



(31) 



(n) Jn+1) <l>k+3{A^-'Sn) <^k+2{A'-'Sn+l] 

e^m+m+l)+i,m ^im+l)k+i,m Hk+2{^^Sn) i/fe+l ( A^^„+i) 



^k+3{^''-'Sn)Hk+liA'Sn+l) - ^k+2i^'-'Sn+l)Hk+2{^'Sn) 
Hk+2iA^Sn)Hk+l{A^Sn+l) 

-f^fc+2(A*~^S'„+i)-f/'fc+i(A*+^S'„) 



Hk+2{A^Sn)Hk+i{A^S n+l) 
Besides, from the identity (fTOj) . with ASn instead of 5"^, we get 

(n+l) (n) Hk{A^+^Sn+l) Hk{A^+^Sn) 

e^m+l)k+l,m Hn.+l)k+l,m Hk+l{ASn+l) Hk+l{ASn) 



2,3, ...,m. (32) 



Hk{A^^^Sn+l)Hk+l{ASn) - Hk{A^+''Sn) Hk+l{ASn+l) 

(33) 



Hk+l{ASn+l)Hk+l{ASn) 
Hk+l{A^Sn)Hk{A^+'Sn+l) 



Hk+1 {ASn+l)Hk+l [ASn] 
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From (fTOl) . we get 



.(n+l) _ (n) _ -^fc+2('5'n+l) Hk+2{Sn) 

■-im+l){k+l),ra-^(m+l)(k+l),^ - H k+l{/^'^+^ S n+l) " i/fc+l(A-+l^„ 



-f^fc+2('S'n+i)-f/'fc+i(A™'+^5'„) — iffc+2(5'„)-fffc+i(A"'+^S'„) 



i/fc+i(A™+iS'„+i)iJ,.+i(A™+i^„) ■ 
Finally, from the identity fll6p . we have, for i = 2, . . . ,m, 

^k+2{^'~'Sn+l) ^k+2{^'''Sn) 



Hk+,{Ar-+^Sn+i)Hk+i{A--+^Sn) 

(34) 



Hk+2{A.Sn) Hk+l{A"^ Sn+l 



-in+1) _ (n) ._ 

'(m+l)k+i,m im+l)k+i,m Hk+l{A^Sn+l) Hk+l{A^S, 



^k+2{A'-^Sn+l)Hk+l{A^Sn) - ^k+2{A'-^Sn)Hk+l{A'Sn+l] 

Hk+i{A^Sn+i)Hk+i{A^Sn) 

Hk+l{A^^^ Sn)Hk+l{A^^^ Sn+l) 



(35) 



i7fc+i(A*5„+i)/7fc+i(A*5„) • 
Then, from the formulae (I33p - (|5S]) . we have 

n/ (n+l) (n) X Hk+i{A^Sn)Hk{A'^'^''Sn+i) ITT Hk+i{A^^^Sn)Hk+i{A^ ^'S'„+i) 

^^^^hm+l)k+^,m ^im+l)k+^,m) " " ( AS„,+i )i/,.+i ( A5„,) ii ( A^^„)i7,+i ( A^^„+i) 

/7fc(A-+i5„H.i)i7fc+i(A-+i5„) 

(36) 



i/fc+i(A5„)iffe+i(A™5„+i) ■ 

Comparing f l30p and 0360 . we obtain the rule (jl]) of the multistep e-algorithm where the lower 
index k is replaced by (m + l)k + 1. 

We can also derive the following formula 

n/ (n+l) (n) \ Hk+2{ASn) Hk+l{A'^ Sn+l) TT- i/fc+i(A*+^5'„)iffc+i(A* ^Sn+l) 

^^2 ^"'+'^'+"'" ^('"+1)^+^''-^ " i7,+i(A™+i^„+i)^fc+i(A-+i5„) i| i/fc+i(A^5„)i/fc+i(A*5„+i) 

_ -f^fc+2(AS'n)/Jfc+i(AS'„+i) ^2^^ 



Hk+i{A^Sn)Hk+i{A-^+^S, 



n+l J 



Comparing flM]) and (jSTj), we obtain the rule (jl]) of the multistep e-algorithm where the lower 
index k is replaced by {m + l)k + 2. 
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Besides, we have, for i = 2, . . . , m. 



.(n+1) in) 



j=i+l 



m i—1 

n( {n+i) _ in) )Y\(e'^"^^^ -e^"^ ) 

V^(m+l)fc+j,m ^(m+l)k+j,mJ 1. 1.^ {m.+l){k+l)+j,m ^ {ni+l){k+l)+j,m^ 
j=i+l j=2 

V^(m+l)(fc+l),m. ^(m+l){fc+l),m/vt(m+l)(fc+l)+l,m ^ {m+l)(fc+l)+l,m/ 



Hk+,{A^Sn+l)Hk+l{A^Sn) \\ Hu+2{AlSn+l)Hk+2{A^S, 
Hk+2{ASn)Hk+l{A"^ Sn+l) Hk+2{A'^Sn)Hk+l{A"^'^^ Sn+1 



/7fc+l(A-+15„,+i)fffc+l(A™+15„) Hk+2{ASn+l)Hk+2{^Sn] 



(3^ 



Comparing fl52]) and fl5Sl) . we obtain the rule (jlj) of the multistep e-algorithm with the lower 
index k replaced by {m + l)k + i + l, for z = 2, 3, . . . , m. Therefore, the multistep e-algorithm has 
been derived from the definition (|6]), and ([8]) of the multistep Shanks' transformation, and we 
have the 

Theorem 4 

The multistep Shanks' transformation defined by and Ij^ can be implemented by the 

recursive rules of the multistep e-algorithm, with the initializations ([5]). 



8 An extended discrete Lotka— Volterra system 

Recently, as explained in Section [H it has been shown that integrable systems are closely related to 
numerical algorithms. On one hand, some numerical algorithms are found to be soliton equations. 
For example, one step of the Qi?-algorithm is equivalent to the time evolution of the finite non- 
periodic Toda lattice [39]. The e-algorithm is nothing but the fully-discrete potential KdV 
equation, and the p-algorithm is considered to be the fully-discrete cylindrical KdV equations 
or the Milne-Thomson equation, see [21l[2S],[nilllIlllS]. On the other hand, integrable systems 
can be used for designing new numerical algorithms. For example, the discrete Lotka- Volterra 
system has applications in numerical algorithms for computing singular values [I9l|20lll0], the 
continuous-time Toda equation leads to a new algorithm for computing the Laplace transform of 
a given analytic function [27], and the discrete relativistic Toda molecule equation leads to a new 
Fade approximation algorithm for formal power series |23] . 

In this section, we will show that there exist a Miura transformation between the multistep 

e-algorithm (jl]) and a discrete integrable system. In fact, if we set fo.^'l^m^ j = — ^i"m; 
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then equation (j4]) is transformed into the extended discrete Lotka-Volterra equation 



m—l m—1 



n "'k'-S^+i n - (n) (n+l) ■ (39) 



2 flW 

i=0 «=0 I "1+1 '"+1 

This equation can be considered as the time discretization, for = —1, of 

dt \ n 



^ ' n^fc-^^+i) = n "fc+in+i+i" n ^k-^-v "^=1,2,..., iv=-i,-2,..., (40) 



which is called the extended Lotka-Volterra equation. This equation was first proposed in [28] 
and it was developed in [TB]. Indeed, with = — 1, (HUl) becomes 



d f"i=r' \ 1 1 



. n ^k-ni^+i = - — - — • (41) 

\ i=0 / 2~ 2~ 

Now, consider n as the discretization of t, and replace the derivative in the left hand side of fl4ip 
by the forward difference A acting on n. The left hand side becomes 

m— 1 m—l 

n in+l) _ TT (n) 

i=0 j=0 

Then, replace Cfc in the first term of the right hand side of fj^Tj) by a[.'^\and, in its second term, 
by 4"+'). We get (jSH]). 

Using the relations f l2T|) and f l23il . we obtain the solution of fl39|l 



a 



a 



(n) i _ -^(m+l)fc+l-^(-m+l)fc+l 

{m+l)k+l-^ ~ ~ in+1) _ (n) ~ ~ pn pn+l 

^(m+l)fc+l,m ^{m+l)k+l,m (m+l)fc+2 (m+l)fc 

(n) _ J- _ (m+l)A;+i-'^^(m+l)fc+j 



that is 



{m+l)k+i-^ (n+l) _ in) pn pn+l ' 

^ {m+l)k+i,m. ^{m+l)k+i,m {m+l)k+i+l {m+l)k+i~l 

z = 2, . . . ,m + 1, 



(n) ^ i7,(A'"+^gOi7,(A-+^5„+0 

(n) _ ( A5'n)i/fc+l ( A5'n+i ) , , 

«(„+l),_:r^+l - Hu+l{^^Sn)Hu{/^rn+lS^^^y ^ > 

(n) ^ Hk+l{A^ Sn)Hk+l{A^ Sn+l) 
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where j = 2, . . . , m, and k = —m + 1, — m + 2, . . ., with the initial values 

2. m-l — OO, , „i_i = ■ ■ ■ = o|_-,''_m-l = m-1 = (45) 

m 2 m-M 2 -'-2 2 

The difference equation fl5^ . with the initial values fHS]) . is said to be the integrable time 
discretization of the extended Lotka-Volterra equation fHU]) in the sense that its solution is given 
by (I42l) - (l4^ . Conversely, the extended discrete Lotka-Volterra equation (HOjl can be seen as the 
time continuation of (l39l) with the initializations fl45l). 



Consider the particular case m = 1. Then, (j4j) reduces to the e-algorithm and equation (139 
becomes 

(n+l) (n) _ 1 1 



(n) (n+l) 



By the dependent variable transformation 



(n+l) (n) 

we obtain the discrete Lotka-Volterra equation 



4""^Mi + <T =<'Mi + -L:i • (46) 



Then, the e-algorithm can be transformed into the discrete Lotka-Volterra equation f H6|) through 
the following Miura transformation 

(n) ^("+1) _ Jn) 



(n+l) (n+2)_ (n+l)' 

"fc-l ^fc-1,1 ^fc-1,1 



Thus, the ^-algorithm can be considered as the discrete Lotka-Volterra equation (H6|) . and more 
generally, the multistep e-algorithm (jlj) is equivalent to the extended discrete Lotka-Volterra 
equation fl5^ . 



9 Conclusion and future researches 

Starting from the recursive rule (jl]) of the multistep e-algorithm, we first obtained, from the 
Hirota's bilinear method, the coupled relations (12T|) - (12^ . Then, applying the Sylvester's identity 
to the determinants Hk{A^Sn) and $fc(A*S'„), we got the formulae (l271) - (l29i) which express the 
quantities e^"^ as ratios of determinants. Thus, we were able to prove that the e^kmS are defined 
as ratios of determinants, and then to derive the recursive rule (jl]) of the multistep e-algorithm, 
with the initializations ([5j), from the determinantal formulae defining the quantities e^"^. It must 
be noticed that, contrarily to the approaches of [15] and [H], we did not make use of the Schweins' 
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determinantal identity, but only of the Sylvester's one. The difficult point was to find to which 
determinants this identity had to be applied. Then, we showed that the multistep e-algorithm 
was related to an extended discrete Lotka-Volterra system. 

When m = 1, the relations ([8]) disappear, and the Hirota's bilinear method leads to a new proof 
that the e-algorithm of Wynn implements the Shanks' sequence transformation and, reciprocally, 
that the quantities computed by this algorithm are expressed by the ratios of Hankel determinants 
defining the Shanks' transformation. 

The approach developed above could possibly be extended to other nonlinear convergence 
acceleration algorithms such as, for example, the g-difference version of the e-algorithm proposed 
in [T3], or its two generalizations given in [S], or the other one presented in [IS], or the general 
e-algorithm of [12] , or the p-algorithm [16] , and the 7-algorithm which generalizes it [3] . Other 
algorithms related to them, such as the qd, the rj, the and the rs-algorithms, and the g- 
decomposition, could also possibly be treated in a similar way (see [11] for their definitions). 
The quantities computed by these algorithms are all defined as ratios of determinants. These 
extensions, as well as extensions to other acceleration algorithms, will be the subject of future 
works. Let us mention that the confiuent form of the multistep ^-algorithm is studied in [10]. It 
leads to a multistep Lotka-Volterra equation. 
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